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INTRODUCTION 


Background 

Previous  calculations  by  another  contractor  using  a  nonstaggered  numeri¬ 
cal  scheme  for  ballistic  flow  with  no  shell  to  barrel  gap  led  to  peculiar 
results  at  the  contact  corner.  In  order  to  proceed  with  the  systematic  devel¬ 
opment  of  a  multi-dimensional,  multi-phase  ballistics  code  it  is  necessary  to 
resolve  this  matter  and  determine  the  role  played  by  the  shell  to  barrel  gap. 
Because  of  JAYCOR's  experience  in  computational  fluid  dynamics,  this  project 
was  initiated  to  address  the  corner  flow. 

Scope  of  the  Work 

In  order  to  produce  results  on  corner  flow  behavior  with  a  minimum  of 
distraction  by  other  ballistics  issues  it  was  decided  to  limit  the  cases 
studied  and  parameters  varied.  The  initial  scope  called  for  runs  for  a  40  mm 
diameter  barrel  with  and  without  a  1%  gap  and  for  a  20  mm  diameter  barrel  with 
no  gap.  A  two-variable  turbulence  model  was  to  be  used  for  the  40  mm  case.  A 
moving  coordinate  system  was  adopted  to  avoid  the  problem  of  an  expanding 
grid,  with  pressure  boundary  conditions  on  the  open  end  to  be  taken  from  the 
earlier,  nonstaggered  grid  calculations.  Accuracy  of  the  calculations  were  to 
be  checked  with  runs  on  a  20  x  20,  40  x  40,  and  80  x  80  grid. 

After  the  initial  results  were  reported  it  was  decided  to  add  additional 
cases:  20  mm  with  a  1%  gap  at  20  x  20  and  40  x  40  resolution.  At  the  same 
time,  it  was  decided  to  explore  other  outlet  flow  conditions  at  the  top  of  the 
gap  to  force  a  net  bypass  flow  around  the  shell. 

The  results  of  these  calculations  were  to  be  analyzed  from  the  point  of 
view  of  phenomena  occurring  in  the  corner  region. 

Summary  of  results 

The  results  are  discussed  in  detail  in  the  next  section  but  they  can  be 
summarized  as  follows. 

(1)  All  of  the  cases  were  successfully  run  by  the  code  except  when  the 
outflow  boundary  conditions  tried  to  drive  the  interior  flow  supersonic  in  the 
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absolute  frame  of  reference  or  when  the  turbulence  model  gave  aphysical  re¬ 
sults. 

(2)  For  the  outlet  boundary  conditions  considered,  the  flow  near  the 
corner  is  considerably  different  with  and  without  a  gap.  The  net  flow  around 
the  projectile  causes  a  flow  splitting  at  the  corner  with  part  of  the  flow 
feeding  the  developing  boundary  layer  and  part  venting  through  the  gap.  The 
flow  in  the  corner  region  tends  to  be  more  complex,  with  a  recirculation  zone, 
counter  current  streams  (in  the  projectile’s  frame  of  reference),  and  a  bound¬ 
ary  layer  all  occurring  within  a  very  small  space. 

(3)  The  gas  temperature  in  the  neighborhood  of  the  projectile  corner  is 
nearly  equal  in  both  cases,  but  the  wall  temperature  at  that  location  can  be 
as  much  as  150°  higher  with  a  gap.  This  difference  is  due  to  the  wall  heating 
in  the  gap  region  that  feeds  the  corner  flow  and  to  the  reduction  of  cooling 
flow  to  the  boundary  layer  caused  by  flow  splitting. 

(4)  The  corner  flow  is  strongly  affected  by  the  prescribed  boundary  con¬ 
ditions.  The  corner  flow  develops  in  about  the  same  length  of  time  as  the 
overall  ballistic  motion  so  the  start  up  transients  are  an  inseparable  part  of 
the  problem.  The  initial  acceleration  of  the  shell  creates  a  low  pressure  re¬ 
gion  at  the  base  which  is  fed  both  by  the  chamber  and  the  gap  -  this  effect  is 
a  major  difference  between  the  gap  and  no  gap  cases.  The  downflow  from  the  gap 
leads  to  a  corner  eddy  that  persists  for  most  of  the  problem  time  and  strongly 
effects  the  pressure  history  in  that  region.  Flow  reversals  in  the  chamber  due 
to  rarefaction  pulses  also  influence  the  developing  flow  and  may  be  caused  by 
inconsistencies  between  the  prescribed  outlet  pressure  and  the  present  calcu¬ 
lation.  The  sudden  application  of  the  gap  outflow  conditions  leads  to  a  rare¬ 
faction  pulse  in  the  gap  that  influences  the  corner  flow  and  the  assumption  of 
a  uniform  outflow  seems  to  contradict  the  counter  current  flow  (in  the  shell’s 
frame  of  reference)  that  tries  to  develop.  These  effects  may  not  affect  the 
qualitative  nature  of  the  flow  but  they  surely  affect  the  quantitative  detail. 

(5)  Turbulence  does  not  seem  to  have  sufficient  time  to  develop  from 
laminar  flow  and  the  use  of  conventional  models  led  to  spurious  results.  Esti¬ 
mates  were  made  that  showed  there  were  not  enough  e-folding  times  to  reach 
transition  based  on  the  growth  of  laminar  disturbances.  Use  of  a  two-variable 
turbulence  model,  based  on  fully  developed  turbulence  arguments,  gave  unphysi- 
cally  large  production  in  the  corner  region  and  caused  the  calculation  to 
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abort.  On  the  other  hand,  calculations  made  with  molecular  viscosity  alone 
were  stable. 


V 
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DISCUSSION  OF  RESULTS 


Standard  resolution,  20  mm,  no  gap 

The  starting  point  for  the  work  done  under  this  project  is  a  20  x  20  cell 
calculation  of  the  flow  in  a  20  mm  diameter  barrel  with  a  shell  that  allows  no 
gap.  Figure  1  shows  the  arrangement  of  mesh  points;  there  is  a  concentration 
of  cells  near  the  wall  to  adequately  resolve  the  developing  boundary  layer  and 
flow  near  the  corner.  Note  that  in  plotting  the  mesh  the  radial  distance  has 
been  expanded  by  a  factor  of  four  so  that  more  detail  can  be  seen;  even  ex¬ 
panded,  the  finest  cells  can  still  be  barely  resolved. 

At  the  initial  instant  the  gas  and  shell  are  at  rest  and  the  gas  pressure 
is  uniform.  The  high  pressure  begins  to  accelerate  the  shell  which  increases 
the  volume  of  the  firing  chamber  and  lowers  the  pressure.  In  all  calculations 
discussed  in  this  report,  the  boundary  pressure  along  the  bottom  was  pre¬ 
scribed  and  the  shell  motion  calculated  from  Newton's  Law.  The  start  up  pro¬ 
cess,  as  well  as  possible  inconsistencies  between  the  previously  calculated 
boundary  pressure  and  the  present  calculation,  causes  a  small  rarefaction  wave 
to  be  generated  and  bounce  off  the  computational  boundaries  a  few  times  before 
diffusing.  The  time  required  for  the  rarefaction  generated  at  the  base  of  the 
shell  to  propagate  to  the  open  boundary  and  back  is  about  0.1  ms.  This  period 
corresponds  to  the  appearance  of  flow  reversal  in  the  corner  at  0.08  ms  and 
0.2  ms.  Figures  2-5  show  the  flow  oscillation  at  times  0.06,  0.08,  0.12,  and 
0.20  ms.  The  vectors  are  scaled  differently  in  each  Figure  and  the  radial 
velocity  component  is  magnified  in  proportion  to  the  expansion  of  the  radial 
coordinate  itself.  No  more  flow  direction  changes  occur  after  this  time  and, 
except  for  scale  effects,  the  corner  flow  is  typified  by  that  shown  in  Fig¬ 
ure  6:  a  small,  steady  radial  flow  toward  the  wall  that  feeds  a  narrow,  grow¬ 
ing  boundary  layer. 

During  the  transient  the  pressure  is  dominated  by  the  overall  chamber 
pressure,  set  by  the  pressure  boundary  conditions,  and  the  slight  gradient 
away  from  the  shell  that  causes  the  gas  to  lag  the  shell.  See  Figure  7  for  a 
typical  pressure  contour  map  at  0.84  ms.  The  early  gas  dynamics  waves  are  also 
planar  and  do  not  readily  show  up  on  these  plots.  At  late  times  the  pressure 
variation  is  still  essentially  axial,  c.f.  Figure  8.  With  the  exception  of  the 
momentary  pertubation  due  to  the  bouncing  gas  dynamics  waves,  the  pressure 
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histories  at  the  centerline  of  the  shell  (Figure  9)  and  at  the  corner  (Figure 
10)  show  the  same  steady  decrease  due  to  expansion. 

The  temperature,  on  the  other  hand,  evolves  considerably  during  the  mo¬ 
tion  of  the  shell.  At  early  times  it  is  nearly  uniform  and  drops  steadily  as 
the  pressure  drops.  See  Figure  11.  As  the  wall  boundary  wall  heating  sets  in 
and  gradients  develop  along  the  wall.  See  Figure  12.  The  temperatures  are 
cooler  near  the  bullet  because  this  is  where  the  inflow  of  gas  from  the  center 
of  the  chamber  feeds  the  boundary  layer.  At  2.29  ms,  a  30°C  temperature  dif¬ 
ferential  has  developed  between  a  point  near  the  corner  and  the  centerline. 
See  Figure  13.  The  difference  at  the  corner  is  much  less.  The  temperature 
nearest  to  the  corner  shows  the  same  general  decrease  evidenced  in  the  pres¬ 
sure  histories.  See  Figure  14.  The  early  oscillations  are  due  to  essentially 
adiabatic  heating  and  cooling  due  to  the  gas  dynamics  waves. 

Standard  resolution,  20  mm,  1%  gap  case 

Next,  the  calculational  domain  was  extended  along  the  axial  direction  to 
allow  the  resolution  of  a  1%  gap.  Figure  15  shows  the  entire  mesh;  the  axial 
and  radial  scales  have  been  further  distorted  so  that  detail  can  be  seen. 
There  are  three  cells  resolving  the  gap  width  and  twenty  along  its  length.  The 
mesh  below  the  shell  is  identical  to  the  one  used  in  the  no  gap  case. 

The  gap  problem  requires  an  additional  boundary  condition:  the  outflow 
condition  at  the  end  of  the  gap  at  the  front  of  the  shell.  This  condition  can 
only  be  determined  by  making  a  complete  calculation  that  extends  far  in  front 
of  the  shell  where  the  ambient  conditions  are  known.  In  order  not  to  detract 
from  the  immediate  objective  of  this  project,  that  of  studying  the  corner 
flow,  it  was  decided  to  use  an  ad  hoc  assumption.  Several  possibilities  were 
discussed  and  tried  including  setting  the  outlet  velocity  to:  10%  more  than 
the  shell  speed;  the  shell  speed;  50%  of  the  shell  speed;  shell  speed  +  90% 
sonic;  90%  sonic;  and  others.  We  now  discuss  one  series  of  runs. 

Since  the  pressure  drop  from  the  back  to  the  front  of  the  shell  is  quite 
large,  it  was  suggested  that  the  flow  might  be  choked,  that  is,  that  a  stand¬ 
ing  shock  would  form  at  the  outlet.  Since,  however,  the  gap  is  quite  narrow, 
10  4  m,  and  viscosity  is  not  negligible  over  the  time  of  the  problem,  vAt/Ar 
~  1,  then  the  outflow  may  be  sonic  only  over  part  of  the  gap.  This  line  of 
argument  led  to  the  proposition  that  the  outflow  be  90%,  say,  of  sonic.  It  was 
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further  proposed  that  the  outflow  be  relative  to  the  shell.  This  means  that 
when  the  shell  reached  a  speed  of  about  100  m/sec,  the  gas  flow  would  be  sonic 
relative  to  the  absolute  frame  of  reference.  Precisely  at  this  time  the  code 
diverged.  Graphics  from  the  last  few  cycles  show  clearly  what  had  taken  place. 

Figure  16  shows  the  velocity  vectors  in  the  gap  region  just  before  the 
code  diverged.  The  outlet  velocity  has  just  gone  slightly  supersonic  in  the 
absolute  frame  of  reference  and  the  third  row  of  velocities  marks  the  upstream 
edge  of  a  shock  that  has  formed.  Figure  17  shows  the  pressure  contours  at  the 
same  instant  and  the  bunching  together  of  the  values  around  an  internal,  maxi¬ 
mum  marks  the  approximate  shock  location.  Since  the  method  is  not  designed  for 
supersonic  flow,  once  it  is  created  a  nonphysical  runaway  thermodynamic  pro¬ 
cess  occurs.  Figure  18  shows  the  temperature  just  ahead  of  the  shock  in  excess 
of  5000  C  and  this  catastrophic  production  of  internal  energy  is  the  process 
that  eventually  causes  the  code  to  fail.  Clearly  the  gas  velocity  cannot  ex¬ 
ceed  sonic  in  the  absolute  frame  of  reference. 

Next  we  tried  an  outlet  condition  that  was  90%  of  sonic  in  the  absolute 
frame  of  reference.  It,  too,  diverged  shortly  after  the  start  of  the  problem 
for  the  same  reason:  an  internal  shock  formed,  despite  the  fact  that  the  out¬ 
let  remained  subsonic.  The  growth  of  the  wall  boundary  layer  is  sufficient  to 
retard  a  significant  part  of  the  flow  and  force  the  central  part  of  the  gap 
flow  to  become  supersonic.  See  Figure  19.  An  additional  process,  just  starting 
to  appear  at  this  time,  is  the  flow  separation  at  the  inlet  to  the  gap  which 
also  constricts  the  flow  and  forces  an  internal  shock  to  form. 

Finally,  the  choice  of  an  outlet  velocity  of  50%  of  sonic  in  the  absolute 
frame  of  reference  worked.  A  maximum  internal  velocity  of  about  90%  sonic 
occurs  0.3  ms  into  the  problem  which  then  reduces  to  nearly  the  outlet  speed 
at  late  times.  The  outlet  conditions  lead  to  a  upstream  propagating  rarefac¬ 
tion  that  accelerates  the  flow  toward  the  outlet.  This  can  be  seen  in  Figures 
20  and  21.  Figures  22  and  23  show  that  the  developing  wall  boundary  layer 
forces  the  flow  toward  the  shell  and  accelerates  it  to  a  value  nearly  twice 
the  exit  speed.  At  this  time  the  chamber  pressure  is  still  greater  than  the 
exit  as  can  be  seen  in  Figure  24.  Around  0.3  ms  the  chamber  pressure  drops  be¬ 
low  that  of  the  exit  causing  the  flow  to  decelerate.  See  Figure  25.  After 
another  0.3  ms  the  gap  flow  settles  down  to  a  balance  between  a  nearly  fully 
developed  boundary  layer  and  a  weak  pressure  gradient.  As  Figure  26  shows  flow 
reversal  occurs  right  at  the  exit  indicating  that  a  uniform  exit  velocity  is 
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probably  not  correct.  The  flow,  however,  seems  to  recover  quickly  so  that  at 
the  gap  it  is  in  the  proper  direction. 

The  maximum  velocity  histories  for  the  above  cases  are  summarized  in  Fig¬ 
ure  27.  Qualitative  agreement  was  gotten  when  the  mesh  was  refined  by  a  factor 
of  two  with  the  maximum  velocity  nearly  going  sonic. 

The  velocity  patterns  in  the  vicinity  of  the  corner  are  shown  in  Figures 
28-31.  The  low  pressure  generated  by  the  initial  acceleration  of  the  shell  at 
first  draws  flow  in  from  the  gap.  As  the  density  of  the  gas  in  the  gap  de¬ 
creases,  the  pressure  falls  and  backflow  finally  ceases  around  0.1  ms.  Further 
flow  from  the  exit  lowers  the  gap  pressure  below  that  in  the  chamber,  drawing 
gas  into  the  gap.  Finally,  the  wall  boundary  layer  develops  and  flow  splitting 
at  the  corner  occurs.  For  the  outlet  conditions  used  in  this  case  the  radial 
flow  along  the  shell  has  been  greatly  enhanced. 

The  pressure  contours  are  similar  to  the  no  gap  case  where  the  pressure 
was  quite  uniform  within  the  chamber  near  the  corner.  See  Figure  32.  However, 
during  the  development  of  the  radial  flow  into  the  boundary  layer  and  the  gap, 
the  corner  and  centerline  pressures  differ  by  as  much  as  300  atmospheres.  See 
figure  33.  There  was  no  discernible  difference  in  the  no  gap  case.  The  center- 
line  curves  from  the  gap  and  no  gap  cases  overlay  one  another  except  for  the 
early  oscillations.  This  agreement  seems  to  indicate  that  the  chamber  pressure 
is  governed  by  the  conditions  at  the  open  boundary  and  that  the  rarefaction 
wave  that  bounced  around  in  the  first  calculation  has  been  filled  in  with  the 
back  flow  from  the  gap. 

The  temperature  field  evolution  is  also  similar  to  the  no  gap  case:  at 
early  time,  uniform  temperatures  primarily  governed  by  expansion  cooling;  at 
late  times,  viscous  heating  of  the  wall  boundary  layer.  The  temperature,  how¬ 
ever,  varies  from  the  centerline  value  as  the  flow  evolves.  See  Figure  34. 
Since  the  wall  flow  is  carried  in  from  the  gas  ahead  of  the  shell,  the  wall 
temperatures  will  depend  upon  the  outlet  conditions,  which  are  probably  not 
totally  correct  in  these  calculations.  The  corner  of  the  shell  and  the  wall 
could  differ  by  more  than  100°  C,  a  variation  not  seen  in  the  no  gap  case.  See 
Figure  35. 
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High  Resolution,  20  mmt  1%  gap  case 

The  previous  calculation  was  repeated  with  every  cell  halved  in  size  in 
both  dimensions  to  produce  a  40  x  40  grid  in  the  chamber  and  a  40  x  7  grid 
within  the  gap. 

The  expanded  flow  in  the  gap  displayed  the  same  characteristics  as  for 
the  low  reslution  case.  The  wall  boundary  layer  reverses  the  flow  one  cell 
below  the  exit  and  the  layer  occupies  2  cells  in  agreement  with  the  earlier 
results.  The  structure  of  the  upward  flow  is  clearer,  with  the  maximum  occur¬ 
ring  away  from  the  shell  as  it  should  since  a  boundary  layer  must  form  there 
too.  See  Figure  36. 

The  pressure  and  temperature  time  histories  at  the  centerline  and  at  the 
corner  compare  favorably  in  Figures  37-40.  There  are  more  variations  in  the 
higher  resolution  run  near  the  corner.  Part  of  the  reason  for  this  is  that  the 
value  is  closer  to  the  corner  where  the  flow  varies  more  rapidly  and  where  it 
may  remain  unsteady. 

The  flow  field  in  the  corner  region  is  shown  in  Figures  41-44.  The  back 
flow  due  to  the  initial  acceleration  is  now  more  accurately  resolved  and  sepa¬ 
rated  flow  from  the  corner  is  clearly  seen.  The  recirculation  eddy  persists 
after  the  backflow  ceases.  Later,  the  wall  boundary  layer  develops  and  the 
flow  along  the  base  splits  at  the  corner  as  occurred  in  the  lower  resolution 
run. 

Comparison  between  the  runs  at  two  different  resolutions  leads  us  to  con¬ 
clude  that  the  20  x  20  calculation  captures  the  quantitaive  essence  of  the 
flow.  The  exact  details  of  the  motion,  especially  at  the  recirculation  zones, 
are  seen  much  more  clearly  in  the  higher  resolution  case.  It  would  probably  be 
appropriate  to  use  these  results  as  boundary  conditions  to  a  very  fine  mesh 
calculation  if  more  detail  were  required. 

Visualization  of  total  flow  field 

The  aspect  ratio  of  the  geometries  in  this  problem  are  so  large  that  it 
is  impossible  to  view  the  flow  over  any  reasonable  extent  and  maintain  the 
proper  geometric  relationships.  We  have  used  a  distorted  presentation  of  the 
data  in  the  corner  region  to  illustrate  the  flow  pattern  there  with  some  suc¬ 
cess.  In  this  section  we  present  the  entire  calculational  domain,  with  even 
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greater  geometric  distortion,  in  order  to  give  the  reader  a  feeling  for  the 
overall  motion  of  the  fluid.  The  technique  employed  is  to  make  every  velocity 
vector  fill  the  cell  it  occupies  so  that  through  a  head-to-tail  connection  the 
viewer's  eye  is  attracted  to  the  streamline  pattern.  The  penalty  for  this  ap¬ 
proach  is  that  very  small  flows  are  over  emphasized,  however,  the  comprehen¬ 
sion  gained  by  the  overview  seams  to  justify  the  attempt. 

The  whole  field  visualization  for  the  no  gap  case  is  shown  in  Figures  45- 
49.  The  flow  reversals,  seen  earlier  in  the  corner  pictures,  occurs  over  the 
whole  field  and  may  be  due  to  an  inconsistency  between  the  shell  motion  and 
the  prescribed  boundary  pressures.  This  method  of  presentation  has  supressed 
the  boundary  layer  by  making  all  velocities  equal;  actually  most  of  the  flow 
is  in  a  thin  wall  layer. 

Figures  50-54  show  the  whole  field  flow  visualization  for  the  1%  gap  case 
wih  standard  resolution.  The  initial  backflow  and  separation  region  can  be 
seen  as  well  as  the  local  rarefaction  that  develops  in  the  gap.  The  recircula¬ 
tion  eddy  builds  up  and  is  shed  from  the  shell  as  the  flow  back  into  the  gap 
is  established.  Finally,  a  splitting  of  the  flow  along  the  shell  base  is 
established. 

Figures  55-58  show  the  whole  field  flow  visualization  for  the  high  reso¬ 
lution  case.  The  evolution  follows  that  of  the  lower  resolution  calculation 
except  that  more  detail  on  the  region  just  behind  the  shell  is  obtained.  The 
visualizations  and  the  earlier  quantitative  comparisons  of  pressures  and  tem¬ 
peratures  leads  us  to  conclude  that  the  20  x  20  is  satisfactorily  detailed. 
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Figure  1.  Mesh  configuration  for  20  x  20  calculation  in  a  20  mm  diameter 

barrel.  The  radial  direction  has  been  expanded  by  a  factor  of  four 
to  allow  more  detail  to  be  seen. 
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figure  2.  Flow  in  the  corner  between  the  shell  base  and  the  wall  at  0.06  ms. 

The  radial  coordinate  and  radial  velocity  component  has  been  great¬ 
ly  enlarged  for  clarity.  At  this  moment  the  flow  relative  to  the 
shell  is  away  from  the  shell. 
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TIME  =  .80000E-04  V  *  63.6 


Figure  3. 


Figure  4 


Flow  in  the  corner  is  described  in  Figure  2.  Time  is  0.08  ms  and  a 
flow  toward  the  shell  in  the  accelerating  frame  of  reference  has 
developed  due  to  gas  dynamics  effects  originating  from  the  start 
up. 


TIME  *  . 12000E-03  V  *  88.8 


.  Flow  in  the  corner  as  described  in  Figure  2.  Time  is  0.12  ms  and 
flow  is  once  again  away  from  the  shell. 
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Figure  5.  Flow  in  the  corner  as  described  in  Figure  2.  Time  is  0.20  ms  as  a 
final,  weak  flow  reversal  has  occurred  due  to  the  second  bounce  of 
the  rarefaction  wave. 
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Figure  6.  Final  configuration  of  corner  flow  (t  =  1.85  ms).  Flow  moves  slow¬ 
ly  from  the  centerline  toward  the  wall  where  it  feeds  a  strong, 
narrow  boundary  layer.  Recall  that  there  is  about  an  order  of 
magnitude  expansion  in  the  radial  coordinate  and  velocity. 
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PRESSURE  CONTOURS  AT  TIMS  •  .34000E-03 

Figure  7.  Pressure  field  at  0.84  ms.  A  simple  axial  gradient  between  the 

shell  and  the  open  boundary  forms.  The  rarefaction  waves  are  not 
apparent  on  any  of  the  pressure  plots  because  they  are  small  and 
without  radial  variation. 


PRESSURE  CONTOURS  AT  TIME  -  .228S8E-32 

Figure  8.  Pressure  distribution  at  late  times  (t  =  2.29  ms).  Except  for  the 
overall  magnitude,  the  pressure  contours  remain  the  same  throughout 
the  bullet  motion. 
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TIME  HISTORY  AT  CENTERLINE  AT  BASE  OF  PROJECTILE 

Figure  9.  Time  history  of  pressure  at  centerline  of  shell.  The  overshooting 
due  to  the  bouncing  gas  dynamics  waves  can  be  seen  at  1  ms  inter¬ 
vals.  Up  to  6  or  7  bounces  can  be  resolved,  but  only  the  first  two 
are  strong  enough  to  change  the  direction  of  the  pressure  variation 
and  consequently  the  flow. 
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TIME  HISTORY  AT  CORNER  QF  PROJECTILE 

Figure  10.  Time  history  of  the  pressure  at  corner  of  shell.  The  behavior  is 
identical  to  that  on  the  centerline  because  the  flow  direction  is 
essentially  axial. 
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TEMPERATURE  CONTOURS  AT  TIME  -  .90000E-04 

Figure  11.  Temperature  contours  in  the  corner  region  at  early  times;  t  =  0.08 
ms.  Contour  lines  are  every  10°C  so  the  temperature  is  everywhere 
within  10°C  of  2940°C. 


TEMPERATURE  CONTOURS  AT  TIME  -  . 14749E-02 


Figure  12. 


Temperature  variations  due  to  wall  heating  are  starting  to  develop 
axially.  The  temperatures  are  cooler  near  the  bullet,  and  more 
like  the  bulk  gas  temperature,  because  this  is  where  the  boundary 
layer  is  fed  by  the  general  radial  flow. 
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Figure  13.  Growing  thermal  layer  in  corner  at  late  time  (t  =  2.29  ms).  The 
corner  itself  is  still  near  the  bulk  gas  temperature  because  of 
the  gas  inflow  along  the  shell. 
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TIME  HISTORY  AT  CORNER  QF  PROJECTILE 

Figure  14.  Time  history  of  temperature  at  the  corner  of  the  shell.  The 
temperature  essentially  traces  the  gas  pressure  including  the 
early  oscillations  due  to  compressive  heating  and  expansion  cool¬ 
ing  caused  by  the  bouncing  gas  dynamics  waves. 
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Figure  15.  Mesh  configuration  for  the  standard  resolution  calculation  of  flow 
with  a  1%  gap.  There  are  20  cells  along  the  gap  and  three  cells 
across  its  width.  The  mesh  below  the  shell  is  identical  to  that 
used  in  the  first  calculation. 


TIKE  =  . 1G38SE-03  V  =  116.8 


0.003923 


0.009373 


Figure  16. 


Velocity  vectors  in  the  gap  region  on  the  last  cycle  before  the 
code  diverged  when  the  outlet  conditions  were  90%  sonic  plus  shell 
speed.  The  radial  coordinate  has  been  expanded  by  a  factor  of 
2000.  The  top  row  of  vectors  is  at  the  outlet  and  have  just 
reached  the  sonic  speed.  Since  this  is  a  prescribed  row,  not  a 
calculated  one,  this  is  not  the  source  of  trouble,  but  the  third 
row  has  gone  slightly  supersonic  and  exceeded  the  range  of  valid¬ 
ity  of  the  method. 
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PRESSURE  CONTOURS  AT  TIKE  •  . ie38SE-03 

Figure  17.  Pressure  contours  for  the  same  case  described  in  Figure  16.  A 
shock,  with  overshoot  ringing,  has  formed  at  the  middle  of  the 
”1.2"  contour  lines  near  the  exit.  The  shock  and  an  attendant 
rarefaction  are  clearly  seen  by  the  bunching  of  the  contour  lines. 


2B00. 

2800. 


TEMPERATURE  CONTOURS  AT  TIME  •  . 1S38SE-03 


Figure  18.  Temperature  contours  for  the  same  case  described  in  Figure  16.  The 
extremely  high  temperature  spike  is  just  upstream  of  the  shock  as 
expected.  The  value  of  the  temperature  is  meaningless  because  the 
method  is  not  intended  to  capture  shocks. 
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TIME  =  . IG046E-03  V  =  113.7 


Figure  19.  Velocity  vectors  in  the  gap  region  just  before  the  code  diverged 

in  the  case  when  the  outlet  velocity  was  90%  of  sonic  speed  in  the 
absolute  frame  of  reference.  The  radial  coordinate  has  been 
expanded  by  a  factor  of  2000.  The  internal  flow  has  been  sped  up 
to  a  supersonic  value  by  the  construction  of  the  wall  boundary 
layer  and  the  backflow  at  the  gap  inlet. 


TIME  *  .6327 IE-04  V  =  39.4 


Figure  20.  Flow  in  gap  when  outlet  velocity  is  50%  sonic.  Radial  coordinate 
expanded  2000  times.  The  rarefaction  wave  from  the  outlet  has 
propagated  about  1200  m/sec  x  0.5  x  10"^  sec  -  .06  m  and  the  flow 
correspondingly  accelerated  to  near  the  exit  velocity. 
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Figure  21.  Flow  in  gap  under  conditions  described  in  Figure  21.  The  expan¬ 
sion  wave  has  traveled  most  of  the  gap  and  the  flow  is  uniform  in 
the  radial  direction  with  the  maximum  value  at  the  exit. 


TIME  *  .  M7B0E-03  V  =  105.2 


Figure  22.  Flow  in  gap  under  conditions  described  in  Figure  20.  Boundary 

layer  development  along  the  wall  has  forced  the  flow  inward  toward 
the  shell  where  it  is  accelerated  to  a  value  that  exceeds  the  exit 
speed. 
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TIME  =  .  241 23E-GJ3  V  =  165.7 


Figure  23.  Flow  in  gap  under  conditions  described  in  Figure  20.  Maximum  flow 
rate  achieved  due  to  constriction  of  wall  boundary  layer.  Maximum 
velocity  nearly  twice  exit  speed. 
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PRESSURE  CONTOURS  AT  TIME  •  .24123E-C53 


Figure  24.  Pressure  contours  at  same  time  shown  in  Figure  23.  Flow  is  still 
being  accelerated  up  the  gap  due  to  a  higher  pressure  in  the 
firing  chamber. 
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PRESSURE  CONTOURS  AT  TIME  »  .32S13E-03 

Figure  25.  Pressure  contours  in  gap  under  conditions  described  in  Figure  20 
at  a  time  just  after  the  maximum  internal  velocity  is  achieved. 

The  wall  boundary  layer  is  still  growing  and  constricting  the  flow 
but  the  chamber  pressure  is  now  lower  than  the  exit  so  the  flow  is 
being  decelerated.  At  this  moment  the  two  trends  are  nearly  in 
balance  so  the  maximum  is  held. 


TIME  «  . 18488E-02  V  ■  518. B 


Figure  26.  Flow  in  the  gap  at  late  times  under  conditions  described  in  Figure 
20.  Boundary  layer  is  only  growing  slowly  and  pressure  drop  is 
nearly  zero  so  that  maximum  flow  velocity  has  settled  on  about  the 
outlet  value.  Note  that  flow  reverses  at  first  cell  from  the 
inlet,  indicating  that  uniform  outlet  velocity  is  probably  not 
correct. 
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31 


Figure  27.  Maximum  velocity  in  the  gap,  not  counting  the  applied  outlet 
value,  for  the  cases  discussed  in  the  previous  figures. 


TIME  *  .  538)7 1E-8H  V  =  33.4 


Figure  28. 


Flow  in  the  vicinity  of  the  corner  for  20  ram  diameter  barrel  with 
1%  gap.  Radial  coordinate  and  velocity  proportionately  exaggerat¬ 
ed.  At  early  times  acceleration  of  shell  causes  a  local  flow 
pressure  to  develop  at  the  base  of  the  shell  which  draws  gas  in 
from  the  gap. 


TIME  =  .943I2E-04  V  =  68.3 


Figure  29.  Flow  in  vicinity  of  corner  for  same  conditions  described  in  Figure 
2  .  After  gas  is  drawn  out  of  gap  the  pressure  lowers  to  a  point 
in  equilibrium  with  the  chamber  and  no  further  flow  occurs. 
Recirculation  eddy  generated  by  initial  motion  persists. 
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TIME  =  • 1 1493E-03  V  «  83.1 


Figure  30.  Flow  in  the  vicinity  of  the  corner  for  the  same  conditions  de¬ 
scribed  in  Figure  28.  As  gas  is  vented  from  exit  of  gap,  pressure 
is  lowered  below  that  in  the  chamber  and  flow  is  pulled  into  the 
gap.  The  radial  flow  across  the  bottom  of  the  shell  splits  into 
part  moving  in  the  wall  boundary  layer  and  part  entering  the  gap. 


TIME  *  .93431E-33  V  *  437.7 


Figure  31.  Flow  in  the  vicinity  of  the  corner  for  the  same  conditions  de¬ 
scribed  in  Figure  28.  The  flow  pattern  is  now  fully  developed. 
It  is  similar  to  that  with  no  gap  except  that  the  radial  flow  is 
increased  to  be  able  to  supply  that  leaking  around  the  shell. 
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PRESSURE  CONTOURS  AT  TIMS  -  . 19488E-02 

Figure  32.  Pressure  contours  in  the  vicinity  of  the  corner  for  the  case 
described  in  Figure  28.  At  late  times  the  pressure  is  nearly 
uniform  within  the  chambers,  with  the  gap  being  at  a  lower  value. 


TIME! SEC) 


TIME  HISTORY  AT  CORNER  OF  PROJECTILE 

Figure  33.  Time  history  of  pressure  at  the  corner  of  shell  and  centerline  for 

20  mm  diameter  barrel  with  1%  gap  (  corner;  -  centerline). 

Because  of  the  expansion  flow  into  the  gap,  the  corner  has  a  lower 
pressure  than  the  centerline  until  equilibrium  is  established. 

This  is  in  contrast  to  the  no  gap  case  where  the  pressures  were 
identical. 
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TIME(SEC) 

TIME  HISTORY  AT  CORNER  OF  PROJECTILE 

Figure  34.  Time  history  of  temperature  at  the  corner  of  the  shell  and  at  the 
centerline  for  a  20  mm  diameter  barrel  with  a  1%  gap.  (  cor¬ 
ner;  - —  centerline).  The  complexity  of  the  flowfield  near  the 

corner  affects  the  viscous  heating  which  varies  considerably  as 
the  flow  develops.  Because  the  wall  flow  originates  from  gas 
ahead  of  the  shell,  it  will  strongly  depend  on  the  outlet  boundary 
conditions . 


Figure  35.  Temperature  contours  at  late  time  for  20  mm  barrel  with  a  1%  gap. 

Viscous  heating  in  the  wall  boundary  layer  dominates  the  field, 
but  the  wall  temperature  at  the  wall  is  150°C  higher  than  the  no 
gap  case  because  the  boundary  layer  is  not  totally  derived  from 
the  circulating  chamber  flow.  The  temperature  at  the  corner  of 
the  shell,  however,  is  nearly  the  same  as  in  the  no  gap  case. 
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TIME  =  .89B23E-03  V  =  403.3 


— t — 

t 

! 

1 

—r~ 

t 

r 

t 

t 

t 

t 

t 

■ 

* 

r 

t 

t 

t 

♦ 

t 

t 

t 

t 

1 

i 

r 

t 

t 

1 

• 

4 

t 

t 

t 

1 

i 

4 

t 

t 

t 

1 

i 

4 

t 

t 

t 

1 

» 

4 

4 

t 

t 

t 

i 

t 

t 

t 

1 

i 

4 

r 

t 

t 

♦ 

t 

4 

t 

t 

1 

1 

i 

4 

t 

t 

t 

1 

• 

4 

t 

t 

t 

I 

* 

4 

t 

t 

t 

* 

* 

4 

♦ 

» 

% 

» 

4 

t 

1 

' 

* 

t 

4 

T 

• 

-  ... 

* 

_ 4 _ 

0.00990 


0.00395 


Figure  36.  High  resolution  calculation  of  flow  in  the  gap  for  a  20  mm  barrel 
with  a  1%  gap.  Boundary  layers  on  both  the  wall  and  the  shell  can 
be  seen. 


TIME(SEC) 

TIME  HISTORY  AT  CORNER  OF  PROJECTILE 

Figure  37.  Time  history  of  pressure  at  corner  of  shell  for  20  mm  barrel  with 

1%  gap  (  40  x  40  mesh;  -  20  x  20  mesh).  Differences  are 

primarily  due  to  details  of  the  flow  about  the  corner  during  the 
backflow  and  reentry  stages  of  development. 
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TIME  HISTORY  AT  CENTERLINE  AT  3ASE  OF  PROJECTILE 
Figure  38.  Time  history  of  pressure  at  base  on  centerline  for  20  mm  barrel 

with  l/o  gap  (  40  x  40  mesh; -  20  x  20  mesh).  Results  are 

essentially  the  same  and  differ  only  in  the  magnitude  of  the  early 
time  gas  dynamics  oscillations. 


TIME(SEC) 

TIME  HISTORY  AT  CORNER  OF  PROJECTILE 

Figure  39.  Time  history  of  temperature  at  the  corner  of  shell  for  a  20  mm 

barrel  with  a  1%  gap  (  40  x  40  mesh;  -  20  x  20  mesh).  As 

with  the  pressure,  the  differences  lie  in  the  details  of  the  flow 
field  which  governs  the  viscous  heating  rate.  The  wall  tempera¬ 
tures  at  the  corner  could  be  150 °C  higher  than  in  the  no  gap  case. 
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TIME  HISTORY  AT  CENTERLINE  AT  BASE  OF  PROJECTILE 

Figure  40.  Time  history  of  temperature  at  base  centerline  for  20  mm  barrel 

with  a  1%  gap  ( -  40  x  40  mesh; -  20  x  20  mesh).  Results  are 

identical  with  the  no  gap  case. 


TIME  >  . 83446E-04  V  *  60.8 


Figure  41.  Flow  pattern  in  the  vicinity  of  the  corner  for  the  case  of  a  20  mm 
barrel  with  a  1%  gap.  The  radial  coordinate  has  been  greatly 
expanded.  Early  time  flow  shows  the  backflow  from  the  gap  into  the 
rarefaction  created  by  the  shell's  acceleration.  In  this  higher 
resolution  calculation  the  separated  flow  from  the  corner  is 
clearly  seen. 
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TIME  =  . 12435E-03  V  *  89.2 


Figure  42.  Flow  pattern  in  the  vicinity  of  the  corner  for  the  case  described 
in  Figure  1.  This  is  the  moment  that  the  pressure  in  the  gap  has 
dropped  to  that  of  the  chamber  and  reentry  flow  is  about  to  be¬ 
gin,  The  recirculation  vortex  caused  by  the  backflow  is  well 
developed. 


TIME  =  • 15120E-03  V  =  107.1 


Figure  43.  Flow  pattern  in  the  vicinity  of  the  corner  for  the  case  described 
in  Figure  41.  Flow  across  the  shell  base  is  established  and  it 
splits  into  the  gap  and  wall  layer  streams. 


39 


TIME  *  .S8737E-03  V  a  348.0 


Figure  44.  Flow  pattern  in  the  vicinity  of  the  corner  for  the  case  described 
in  Figure  41.  The  wall  boundary  layer  is  established  and  well 
resolved. 


TIKE  »  . 60000E-04  V  a  47.9 


Figure  45.  Whole  field  visualization  of  flow  in  20  mm  no  gap  case.  All 

vectors  have  been  expanded  to  fill  their  computational  cell  and 
all  cells  have  been  expanded  to  a  square  array.  At  this  early 
instant  flow  is  leaving  the  computational  region  as  the  gas  lags 
the  shell’s  motion. 
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TIME  =  . 80000E-04  V  =  83.6 


Figure  46.  See  caption  for  Figure  45.  Flow  reversal  occurs  throughout  the 
field,  probably  the  result  of  a  rarefaction  wave  set  up  at  the 
shell.  An  inconsistency  between  the  calculated  shell  velocity  and 
the  prescribed  pressure  boundary  condition  may  be  a  source. 


TIME  s  . I2G30E-03  V  =  88.8 


Figure  47.  See  caption  for  figure  45.  Flow  returns  to  normal  direction  as 
rarefaction  wave  leaves  system. 
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TIME  *  .20000E-03  V  3  MG. 3 


Figure  48.  See  caption  for  Figure  45.  A  second  flow  reversal  accompanies 
another  rarefaction  wave. 


TIME  *  . 22Q9SE-02  V  =  648.0 


Figure  49.  See  caption  for  Figure  45.  At  late  times  flow  pattern  settles  down 
to  an  outflow  expansion  which  is  directed  slightly  toward  the 
outer  wall  to  feed  the  wall  boundary  layer. 
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TIME  =  .53071E-04  V  =  39.4 


Figure  50.  Whole  field  visualization  for  20  mm,  1%  gap  case  with  standard 

resolution.  The  backflow  from  the  gap  can  be  seen  and  the  first 
stages  of  a  separation  region.  A  local  rarefaction  in  the  gap 
develops  as  flow  exits  both  ends  of  the  gap. 


TIME  =  .94312E-04  V  =  68.6 


Figure  51.  See  caption  for  Figure  50.  The  flow  separation  region  at  the 
corner  is  well  developed  and  the  wall  boundary  layer  has  been 
established  up  to  the  corner 
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Figure  52.  See  caption  for  Figure  50.  The  recirculation  eddy  has  been  pushed 
from  the  base  of  the  shell  as  flow  into  the  gap  is  established.  A 
local  rarefaction  has  formed  in  the  fifth  row  below  the  shell. 


TIME  *  , 18846E-03  V  =  131.9 


Figure  53.  See  caption  for  Figure  50.  Downward  flow  is  being  reestablished 

by  the  upper  propagating  rarefaction  wave.  The  separation  eddy  is 
smaller  and  nearly  diffused. 
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TIKE  =  . 11668E-02  V  =  451 ,0 


Figure  54,  See  caption  for  Figure  50.  Final  flow  pattern  is  established  with 
flow  splitting  of  the  base  flow  into  the  gap  and  the  wall  boundary 
layer. 


TIME  =  . 6G861E-04  V  *  48.9 


Figure  55.  Whole  field  visualization  for  the  20  mm,  10%  gap  case  with  high 
resolution.  The  backflow  from  the  gap  has  begun  with  a  small 
recirculation  just  forming  along  the  base.  A  rarefaction  has 
formed  in  the  gap  due  to  the  outflow  at  both  ends. 
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TIME  =  .13412E-33  V  «  9B.8 


Figure  56.  See  caption  for  Figure  55.  The  recirculation  region  at  the  base 

of  the  shell  has  expanded  and  flow  splitting  at  the  gap  has  begun. 


TIME  =  .  18710E-C3  V  =  13(3.3 


Figure  57.  See  caption  for  Figure  55.  Flow  reversal  in  the  chamber  has 

occurred  and  the  recirculation  eddy  has  been  shed  from  the  base  of 
the  shell. 
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TIME  =  .89623E-03  V  *  403.3 


Figure  58.  See  caption  for  Figure  55.  The  flow  is  approaching  its  final 

pattern.  The  expansion  flow  is  away  from  the  shell  and  the  flow 
splitting  into  the  gap  and  the  boundary  layer  are  clearly  seen. 
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MATHEMATICAL  MODELS  AND  NUMERICAL  ALGORITHMS 
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NOMENCLATURE 


2 

a  acceleration  of  the  projectile  (m/sec  ) 

E  energy  (J/kg) 

I  index  of  rightmost  column 

i  index  denoting  cell  columns 

J  index  of  topmost  row 

j  index  denoting  cell  rows 

M  mass  of  the  projectile  (kg) 

n  index  denoting  time  level 

p  pressure  (N/m^) 

r  radial  distance  (m) 

T  temperature 

t  time  (sec) 

5t  time  increment  (sec) 

u  radial  velocity  (m/sec) 

v  axial  velocity  (m/sec) 

V  projectile  velocity  (m/sec) 

z  axial  distance  (m) 

p  density  (kg/m  ) 

p  viscosity  of  fluid 

X  viscosity  coefficient 

y  ratio  of  specific  heats,  value  used  was  1.27118. 

o 

a  turbulent  thermal  diffusivity  (m  /sec) 

o 

AA.  •  area  of  computational  cell  i,j  (m  ) 

1  >  J 
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GOVERNING  EQUATIONS 


The  equations  assumed  to  describe  the  expanding  flow  within  the  barrel  of 
the  gun  are  the  compressible  Navier-Stokes  equations  written  in  axisymmetric 
cylindrical  polar  coordinates,  plus  an  appropriate  equation  of  state  to  close 
the  system.  These  equations  are: 

Continuity: 


lf  +  7  37  (pur)  +h  (pv)  =  0 


(i) 


Radial  momentum: 


_3_ 

3t 


(  pu)  + 


r  9r  *• 


pu2r)+|_(Puv)  =  -|£  +  14r(yr  £)-lHu 


r  3r 


3r 


(2) 


Axial  momentum: 


If  (pv)  +7  77  (puvr)  +h  (pv 2 ) =  'If +  7If  M 


9u  ,  3v 
8z  ITr 


)] 


+  2 


3z 


3_ 

3z 


(3) 


Energy: 


If  (oE)  +  f!f  (rupE)  +TF 


<’-) ' I E? If  (™ + 


2  2 
,  3u  ,  3v  3u 

+  y  2  37  +-37  +  -3F 


3  2  2 

+  2  p-  +  2  - 

3z  r 


A  A 
(4) 
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where 


(5) 


Nobel-Abel  Equation  of  State: 


p  =  p(p,E) 


ec  =  0.0010775  m3/kg 


(6) 


Ballistic  Calculation: 


F  =  Ma 


(7) 


The  viscosity  coefficient,  y,  is  either  calculated  from  Sutherland  formula 


4.0608  x  10  8  E3/2 

E  +  1.4232  x  105 


y 


or  y  from  the  k-e  turbulence  model  of  Pope  and  Whitelaw,  or  y  is  taken  as  a 
constant.  The  viscosity  coefficient  X  is  equal  to  -2/3  y. 

Due  to  the  desire  to  keep  the  calculation  domain  constant,  a  transforma¬ 
tion  to  the  coordinate  system  fixed  with  respect  to  the  moving  shell  is  incor¬ 
porated  into  these  equations  in  the  following  manner.  In  the  barrel  fixed  sys¬ 
tem,  the  shell  moves  past  stationary  walls  with  a  velocity  V,  creating  a  gas 
velocity  v.  This  is  the  axial  velocity  in  Eqs.  (l)-(5).  In  the  shell  fixed 
frame  of  reference,  the  shell  is  stationary,  but  the  gun  barrel  now  moves  past 
the  shell  with  velocity  -V.  The  axial  gas  velocity  seen  in  this  frame  is  vfe  = 
v  -  V.  Hence,  the  transformation  of  the  equation  can  be  made  by 


t 


z^  =  z  -  /  V(t")  dt 


0 


vb  -  v  -  V 
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If  this  transformation  is  used  in  Eqs.  (l)-(5),  the  only  change  which  occurs 
is  seen  in  the  axial  momentum  equation;  all  the  other  equations  remain  the 
same,  but  with  the  subscript  b  appended  to  each  variable.  The  axial  momentum 
equation  is  altered  by  the  addition  of  one  term  which  accounts  for  the  un¬ 
steady  motion  of  the  shell.  Dropping  all  the  subscripts  which  appear  in  the 
equation  due  to  the  transformation,  the  axial  momentum  equation  in  the  shell 
fixed  coordinate  system  is 


_3_ 

3t 


(pv>  +  7  I7  (puvr>  +  -57  (pv2)  = 


3z 


“*  P 


dV  3p  ,13  r  r  3u  ,  3  v  ^  i 


dt 


^3z  3r- 


+ 


*  i  (>■£)  + -Ir 


3z 


3z 


(3-A) 


Equations  (l)-(7),  with  Eq.  (3)  replaced  by  Eq.  (3-A),  are  the  governing  equa¬ 
tions  solved  by  JAYCOR's  SAGITA  code. 
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FINITE  DIFFERENCE  EQUATIONS 


An  extended  version  of  the  Implicit  Continuous-fluid  Eulerian  (ICE)  tech¬ 
nique  of  Harlow  and  Amsden^  for  compressible  flow  is  used  in  JAYCOR's  SAGITA 

(2) 

code  which  is  one  application  of  the  EITACC  methodology. 

A  fragment  of  the  Eulerian  finite-difference  mesh  is  shown  in  Fig.  1, 
which  illustrates  the  centering  of  the  field  variables  relative  to  a  typical 
cell.  The  index  i  counts  cell  centers  from  left  to  right,  while  j  counts  them 
from  bottom  to  top.  6r^  and  are  the  lengths  of  a  cell. 


Figure  1.  Layout  of  Variables  and  Indices 
in  the  Mesh. 

The  form  of  the  finite  difference  equations  are  as  follows: 

Ballistic  Calculation: 


rn+l  ,.n  .  6t  r  n+1  .. 
/  =  V  +  —  )  p.  .  AA 

M  .  .  i,i 

1.3 


i.3 


(8) 


Continuity: 


n+1 
P.  , 


Pn  .  +  fit 
1.3 


3p 


n+1 


at  . 

1,3 


(9) 
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Momentum: 


(pUCj  =  ^PU^i , j  +  6t  if  (pu) 


i,  3 


_  i£ 

3r 


n+1 


i»3 


(10) 


=  (pv)i,j  +  6t  If  (pv) 


3_ 

3t 


1  >  3 


_  i£ 

3z 


n+1 


i,3 


(ID 


Energy : 


(pE)±tj =  ^pE^,j  +  6t  If  (pE)  .  .  (12> 

1  >  J 

where  the  index  n+1  indicates  the  advanced-time  value.  The  [  ]  means  that  the 
pressure  is  not  included  in  this  momentum  derivative  terra. 


In  order  to  eliminate  the  sonic  speed  stability  limitation,  the  pressure 
must  be  implicitly  included  in  the  calculation  procedure.  This  is  accomplished 
by  incorporating  the  values  of  the  flux  at  the  forward  time  step,  Eqs.  (10) 
and  (11),  into  the  continuity  equation  (1)  as  follows;  In  shorthand  vector 
notation,  where  w  =  (u,v),  we  have 


_3p 

3t 


+  div(pw)?+^  =  0 
i,3  1,2 


or 


TT  +  div[pw]”+*  -  6t  div  grad  p?+*.  =  ( 
dt  i  j  1»3  G  *i ,3 


(13) 


This  gives  an  equation  for  the  new  pressures,  p?+^ ,  when  (10)  and  (11)  are 

1  )  J 

substituted  in  (13). 
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SOLUTION  PROCEDURE 


The  specific  steps  which  comprise  one  calculation  cycle  are  as  follows: 

1.  From  known  values  at  time  step  n,  and  the  latest  iterated  values  at 
n+1  compute 

* 

Vn+1  from  Eq.  (8) 

from  Eq.  (9) 

without  the  5t  -7^7  term  from  Eq.  (10) 
i»  J 

n+1 

term  from  Eq.  (11) 

i,  j 

(pE)n+*.  from  Eq.  (12) 

2.  These  values  at  n+1  are  placed  in  the  Poisson  Eq.  (13).  This  equation 

is  then  solved  for  p?+^  by  a  direct  method,  rather  than  SOR,  to  avoid  slow 

or  non-convergence  problems  caused  by  high  aspect  ratio  finite  difference 
meshes • 

3.  These  values  of  p?+^  are  used  to  complete  the  calculations  of  (pu)?~^ 

+1  1 5  J  x » J 

and  (pv)?  .  in  Eqs.  (10)  and  (11). 

1  >  1 

4.  After  steps  1-3,  V,  p,  pu,  pv,  pE,  p  are  now  defined  at  the  new  time 
step  n+1.  These  are  compared  with  the  previous  iterated  values.  If  convergence 
is  obtained  then  continue  to  the  next  calculation  cycle;  otherwise,  repeat  the 
above  steps. 


r  ^n+1 

'  PV  i> j 


without  the  6t 

3z 


n+1 

Pi>  j 


r  ^n+1 
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INITIAL  CONDITIONS 


The  initial  conditions  were  taken  from  the  BRL  computer  run  output  sent 
as  an  attachment  to  ARRADCOM  letter  dated  20  December  1979.  These  initial  con¬ 
ditions  are  as  follows: 

8  2 

p  =  3.0  x  10°  N/m  everywhere  in  the  computational  mesh. 

3 

p  -  218.66  kg/m  everywhere  in  the  computational  mesh. 

pE  in  J/m  is  computed  from  the  equation  of  state  for  every  cell  in  the 
computational  mesh. 

O 

pu  =  0.0  kg/m  -sec  everywhere  in  the  computational  mesh, 
pv  =  0.0  everywhere  in  the  computational  mesh. 
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BOUNDARY  CONDITIONS 


Boundary  conditions  are  derived  in  such  a  way  as  to  describe  the  various 
walls  or  surfaces  that  may  occur.  Below  is  the  description  of  the  various 
boundary  conditions  used  at  the  centerline  (left  of  the  computational  mesh), 
at  the  wall  of  the  gun  (right  of  the  computational  mesh),  4  cm  away  from  the 
projectile  base  (bottom  of  computational  mesh),  and  at  the  projectile  base 
(top  of  computational  mesh).  I  and  J  represent  the  total  number  of  columns  and 
rows  used,  with  columns  i  =  1  and  I  and  j  =  1  and  J  being  boundary  cells  out¬ 
side  the  physical  domain. 


Boundary  Conditions  on  Centerline: 


pl,j  '  P2,j-  »1,J  ’  »2,j>  -  (PE)2,J. 

(pu)2,j  =  0,  (Pv)1}j  =  (Pv)2>j 


Boundary  Conditions  on  the  Wall  of  Gun: 

pi,j  =  Pi-i,j>  pi,j  =  pi-i,j’  (pE)i,j  =  <pE>i-i,j’ 

(  Pu ) x  ^  j  =  0, 

(pv)T  •  =  (pv)  where  (pv)  is  defined  such  that  the  wall  is  moving  at 

1  9  J  W  W 

velocity  -  V, 


Boundary  Conditions  4  cm  from  the  Projectile  Base: 


Pi  1  =  Pb  where  Pb  is  taken  from  the  BRL  computer  run. 


Pj^  j,  ( pE )^  j,  (Pu)i}1,  (pv)isj  are  all  linearly  extrapolated. 
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Boundary  Conditions  at  the  Projectile  Base: 


pi,J  "  pi,J  =  pi,J-l>  <pE>i,J  =  (PE)i,J-l> 

(pu^i,J  =  ~  ^pu^i, J-l >  (pv)i,J  =  0 

Boundary  Conditions  next  to  bullet  for  Gap: 


T.j  ‘  Pb,j' 

pb,j  ■  pb,j-  <«*%,!  ■  <pE>b,j’ 

(  pu  )—  =  0 . , 

D  >  J 

(pv)¥jJ  -  -  (pv)bJ 

Boundary  Conditions 

for  top  of  Gap: 

pi,J  =  pi , J-l » 

pi,J  =  pi,J-l>  (pE)i,J  =  (PE)i,J-l» 

(pu)i}J  =  ( Pu)lj j-i 
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TURBULENCE  MODEL 


The  model  I  two-equation  k-e  turbulence  model  as  used  by  Pope  and  White- 
law^)  was  implemented  in  the  code  as  follows: 

9k  .  9k  9k  19  ,  9k  ,  9  3k  ,  D 

-T—  +U-T-+V-T—  =  —  rv^—  +-T—  v-r—  +  P  -  e 

9t  9r  9z  r  9r  9r  9z  9z 


+  V  ft*  +  P  (CHP  Cf?e) 


3e  3c  ,  3e  _  1  3  3e 

3t  U  3r  V  3z  r  3r  r  3r  3z  3z  k  ^elA 


e2 


2 

v  =  minfc  k  /e,  v  n  +  c  v.fr  -  r)l 
L  p  mol  w  * v  w  ' J 


P 


v 


2  2 
9u  ,  jlv  du_ 

9z  9r  9r 


u_ 

r 


2 


c  =  0.09 
V 

c  ,  =  1.45 
el 


c 


e2 


1.90 


c  =  0.04 
w 


V 

* 


j 


This  model  was  used  to  run  the  problem  under  task  F.3.1.5,  which  is  the 
40  mm  gap  problem.  This  model  gave  a  viscosity,  v,  about  2  or  3  orders  of  mag¬ 
nitude  larger  than  vmol  +  cwv*(rw  -  r)  where  vmol  is  the  viscosity  from  Suth¬ 
erland  equation. 
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APPENDIX  B 


DISCUSSION  OF  ADDITIONAL  CALCULATIONS 
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Below  is  a  summary  of  the  calculations  called  out  in  the  contract  (under¬ 
lined)  and  additional  calculations  performed  by  JAYCOR  to  clarify  the  behav¬ 
ior.  These  runs  were  made  with  no  net  bypass  flow,  which  was  later  decided  not 
to  be  the  preferred  boundary  conditions  thereby  necessitating  the  runs  de¬ 
scribed  in  the  main  body  of  the  text. 

F.  3. 1 . 1 .  No  gap  problem  with  20  mm  gun.  The  computational  grid  spacing 
used  was  the  same  as  used  by  the  BRL  computer  run  dated  20  December  1979  in 
the  radial  direction,  and  a  fixed  increment  of  0.002  in  axial  direction. 

F.3. 1.2.  No  gap  problem  with  40  mm  gun.  Two  (2)  runs  were  made  for  this 
task  as  follows: 

F.  3. 1.2.  A.  No  gap  problem  with  40  mm  gun  and  projectile  weight  of  480 
grams.  This  weight  was  used  because  it  scales  with  the  20  mm  gun  with  weight 
of  120  grams.  This  was  done  so  comparisons  could  be  made  to  the  BRL  computer 
run  dated  20  December  1979.  The  computational  grid  spacing  used  was  the  same 
as  m  F.3.1.1. 

F.3.I.2.B.  No  gap  problem  with  40  mm  gun  and  projectile  weight  of  960 
grams.  The  computational  grid  spacing  used  was  the  same  as  in  F.3.1.1. 

F.3.1.3.  No  gap  problem  with  40  mm  gun  and  turbulence  model.  The  projec¬ 
tile  weight  used  was  480  grams  so  again  comparisons  could  be  made  to  the  BRL 
computer  run.  This  was  run  with  a  fixed  viscosity  of  10  ^  [about  two  (2) 
orders  of  magnitude  larger  than  from  Sutherland  equation).  The  computational 
grid  spacing  used  was  the  same  as  in  F.3.1.1. 

F.3. 1.4.  Gap  problem  with  40  mm  gun.  The  computational  grid  spacing  was 
the  same  as  the  BRL  computer  run  in  the  radial  direction,  and  using  an  incre¬ 
ment  of  0.002  meters  up  to  three  (3)  cells  above  the  gap  and  then  an  increment 
of  0.02  meters  to  the  top  of  the  gap.  Three  (3)  runs  were  made  for  this  task 
as  follows: 

F.  3. 1.4. A.  Gap  problem  with  20  mm  gun.  This  run  used  as  the  boundary 
condition  for  (pv)  at  the  top  of  the  gap  to  vary  linearly  from  the  projectile 
(v  =  0)  to  the  gun  wall  (v  =  -V)  as  discussed  in  the  boundary  condition 
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section.  Note  that  this  run  hit  time  limit  at  t  =  0.0019305  and  on  the  graphs 
this  time  is  plotted  with  the  BRL  computer  run  at  t  =  0.0022917. 

F.3.I.4.B.  Gap  problem  with  40  mm  gun  and  projectile  weight  of  480 

grams.  The  weight  was  changed  so  to  scale  with  the  BRL  computer  run.  This  run 
used  as  the  boundary  condition  for  (pv)  at  the  top  of  the  gap  to  correspond  to 
the  velocity  of  the  projectile. 

F.3.I.4.C.  Gap  problem  with  40  ram  gun  and  projectile  weight  of  960 

grams.  This  run  used  the  same  boundary  conditions  for  (pv)  as  F.3.I.4.B. 

F. 3. 1 . 5.  Gap  problem  with  40  mm  gun  and  turbulence  model.  The  computa¬ 
tional  grid  spacing  used  was  the  same  as  in  F.3.1.4.  Three  (3)  runs  were  made 
for  this  task  as  follows: 

F.  3.  1.5. A.  Gap  problem  with  40  mm  gun  with  the  two-equation  k-e  turbu¬ 
lence  model  as  used  by  Pope  and  Whitelaw  and  projectile  weight  of  480  grams. 

—3 

F.3.I.5.B.  Gap  problem  with  40  mm  gun  with  a  fixed  viscosity  of  10 

(about  3  orders  of  magnitude  larger  than  from  Sutherland  equation)  and  projec¬ 
tile  weight  of  480  grams. 

F.3.I.5.C.  Gap  problem  with  40  mm  gun  with  a  fixed  viscosity  of  10”^ 

(about  2  orders  of  magnitude  larger  than  from  Sutherland  equation)  and  projec¬ 
tile  weight  of  480  grams. 

F. 3. 1 . 6.  Gap  problem  with  40  mm  gun,  turbulence  model  and  mesh  refine¬ 
ment.  Two  (2)  runs  were  made  for  this  task  as  follows: 

F. 3. 1.6. A.  Gap  problem  with  40  mm  gun,  fixed  viscosity  of  10”^,  projec¬ 
tile  weight  of  480  grams,  one  grid  halving  in  both  the  radial  and  axial  direc¬ 
tion.  This  run  required  one  (1)  hour  of  CP  time  on  the  CDC  7600  computer. 

F.3.I.6.B.  Gap  problem  with  40  mm  gun,  fixed  viscosity  of  10~\  projec¬ 
tile  weight  of  480  grams,  two  (2)  grid  halvings  in  the  axial  direction  because 
of  the  CP  time  that  would  have  been  required. 

Of  the  twelve  runs  described  above,  some  are  of  greater  interest  in  the 
light  of  the  net  bypass  flow  cases  run  later.  A  discussion  of  the  comparison 
is  given  below  in  order  to  shed  light  on  the  difference  caused  by  changing  the 
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outflow  boundary  conditions  and  to  document  the  experience  gained  with  a  two- 
variable  turbulence  model  and  with  varying  the  radial  resolution. 

The  initial  calculations  were  made  with  zero  net  bypass  flow  around  the 
shell.  Figures  B-l  through  B-5  show  the  whole  field  visualization  of  the  flow 
evolution  in  this  case  which  is  typical  of  the  variation  seen  in  all  cases 
with  no  net  bypass.  The  early  moments  show  the  expansion  region  pulling  gas  in 
from  the  gap.  Recall  that  there  is  considerable  distortion  of  the  cells  and 
velocity  vectors  in  this  representation  and  that  only  the  pattern  is  qualita¬ 
tively  correct.  A  recirculation  develops  but  cannot  be  maintained  although 
reinforced  by  the  flow  reversal  caused  by  the  rarefaction  waves;  this  is  qual¬ 
itatively  the  same  as  in  the  net  bypass  cases.  Figures  50-54.  The  final  flow 
patterns,  Figure  B-5,  are  similar  to  the  bypass  case  in  the  chamber  but  the 
gap  shows  only  downward  flow  in  the  shell’s  frame  of  reference.  The  pressure 
history,  Figure  B-6,  is  the  same  as  before,  but  the  temperature  field,  Figures 
B-7  and  B-8,  is  considerably  cooler  than  with  bypass  due  to  the  much  smaller 
velocity  gradients  in  the  gap.  The  remaining  runs  were  concentrated  on  the  40 
mm  diameter  barrel. 

Figures  B-9  through  B-13  present  results  on  the  40  mm,  no  gap  case  where 
the  shell’s  mass  has  been  scaled  with  the  square  of  its  base  area  to  produce 
the  same  acceleration  as  in  the  20  mm  case.  The  results  for  flow  pattern, 
pressure,  and  temperature  are  essentially  the  same  as  in  the  20  mm  case. 

Figures  B-14  through  B-18  show  the  flow  pattern  in  the  no  gap  case  when 

/  0 

the  fluid  viscosity  is  increased  to  10  ctn  /sec,  about  two  orders  of  magni¬ 
tude  larger  than  the  value  given  by  the  Sutherland  equation.  The  variations  of 
pressure  and  temperature  are  nearly  identical  to  the  lower  viscosity  runs.  The 
high  viscosity  values  were  explored  as  a  substitute  for  a  turbulence  model 
after  problems  developed  with  standard  models. 

Figures  B-l 9  through  B-26  show  the  results  for  the  40  mm  diameter  barrel 
with  a  10%  gap  and  no  net  bypass  flow.  The  same  sequence  of  patterns  that 
occurred  in  the  20  mm  case  reappears:  initial  expansion  and  circulation  in  the 
wake  of  the  shell  dissipation  and  migration  of  the  eddy;  reestablishment  of 
the  circulation  pattern  due  to  a  rarefaction  wave;  and  the  final,  expansion 
flow  out  the  bottom  boundary.  One  significant  difference  is  that  the  wider  gap 
permits  a  boundary  layer  to  develop  on  the  shell  which  leads  to  countercurrent 
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flow  in  the  gap,  Figure  B-22.  Some  flow  splitting  at  the  shell  corner  does 
occur,  but  it  is  not  the  dominant  feature  as  in  the  net  bypass  flow  situ- 
ations.  The  pressure  variation  is  identical  to  the  20  mra  case  because  the  pre¬ 
scribed  boundary  conditions  dominate  the  field,  but  the  wall  temperatures  are 
slightly  lower  because  the  scale  of  the  chamber  has  been  increased.  Internal, 
thermal  eddys  originating  from  the  corner  flow  can  be  seen  in  the  temperature 
contours,  Figure  B-26. 

When  the  turbulence  model  of  Pope  and  Whitelaw  was  used  the  calculations 
invariably  diverged  after  the  shell  reached  sizeable  speeds.  The  problem 
stemmed  from  a  catastrophic  production  of  turbulence  energy  in  the  neighbor¬ 
hood  of  the  corner.  This,  in  turn,  arose  from  the  high  shear  field  produced  at 
the  separation  point  and  the  strong  rarefaction  effects  of  the  accelerating 
shell.  The  model,  based  on  incompressible  and  parallel  flow  data,  is  evidently 
not  applicable  in  this  circumstance.  It  is  unlikely  that  there  is  a  viable 
theory  for  this  flow,  however  coupled  with  the  short  duration  of  this  problem 
fully  developed  turbulence  probably  does  not  occur.  In  any  event,  to  study  the 
effects  of  "turbulence-like"  diffusion,  we  overwrote  the  calculated  turbulence 
quantities,  in  particular  the  eddy  viscosity,  whenever  they  exceeded  the  val¬ 
ues  encountered  in  pipe  flow  at  the  same  mean  velocity,  pipe  diameter,  and 
distance  from  the  wall.  This  modification  did  allow  the  code  to  run  and  the 
results  are  shown  in  Figure  B-27  through  B-34.  The  flow  fields  are  qualita¬ 
tively  similar,,  although  the  boundary  layer  is  thicker  and  countercurrent  flow 
in  the  gap  has  greater  difficulty  forming.  The  pressure  history  is  again  the 
same,  but  the  corner  temperature  is  higher  by  about  150 °C  due  to  accentuated 
wall  friction. 

—3  2 

An  even  greater  viscosity  of  10  cm  /sec  was  chosen  for  the  calculation 
presented  in  Figures  B-35  through  B-42.  Again,  the  flow  pattern  and  pressure 
is  essentially  the  same  but  wall  heating  has  raised  the  temperature  400 °C 
above  the  normal  laminar  case.  This  variation  in  viscosity  over  almost  three 
orders  of  magnitude  is  fairly  convincing  evidence  that  the  flow  pattern  will 
exist  independent  of  the  exact  form  of  the  diffusion. 

On  the  other  hand,  there  is  a  marked  difference  in  the  corner  flow  with 
and  without  net  bypass  flow.  With  bypass  there  is  a  strong  flow  splitting  at 
the  corner  in  the  20  mm  case  that  is  not  seen  when  the  gas  moves  with  the 
shell.  The  40  mm  barrel  with  gap  does  show  the  effect. 
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Finally,  calculations  were  made  with  20,  40,  and  80  cells  in  the  radial 
direction.  The  results  shown  in  whole  field  visualization  are  found  in  Figures 
B-43  through  B-71.  The  calculations  agree  very  closely  over  the  whole  range 
and  at  all  times,  indicating  that  the  lower  resolution  runs  are  sufficiently 
accurate  for  the  present  purposes.  The  largest  difference  was  a  20°C  change  in 
wall  temperature  between  the  20  x  20  and  the  20  x  80  runs. 
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TIME  =  . 58456E -04  V  =  46.7 


Figure  B-l.  Whole  field  visualization  of  20  mm,  1%  gap  when  no  net  flow 
boundary  conditions  are  used  at  the  exit  of  the  gap.  All  vectors  have 
been  expanded  to  fill  their  computational  cell  and  all  cells  have  been 
expanded  to  a  square  array.  Considerable  geometric  distortion  is  intro¬ 
duced  but  the  flow  directions  are  clearly  seen.  At  this  early  time  a  rare¬ 
faction  has  developed  behind  the  shell  which  is  pulling  in  gas  from  the 
gap.  The  flow  pattern  is  similar  to  that  of  Figure  50  where  net  flow  past 
the  shell  was  imposed, 

TIME  -  . I I6I7E-03  V  =  86.1 


Figure  B-2,  See  caption  for  Figure  B-l.  Flow  separation  eddy  has  grown  in 

the  shell fs  wake  and  is  still  qualitatively  similar  to  the  net  flow  calcu¬ 
lation  shown  in  Figure  51. 
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Figure  B-3.  See  caption  for  Figure  B-l.  The  recirculation  region  has  been 

dispersed  for  the  moment  as  a  local  rarefaction  develops  in  its  place.  By 
this  time  in  the  net  flow  calculation,  flow  up  the  gap  has  developed.  See 
Figure  52. 
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Figure  B-4.  See  caption  for  Figure  B-l.  The  recirculation  has  been  reestab¬ 
lished  at  the  corner  of  the  shell  aided  by  the  flow  reversal  of  the  rare¬ 
faction  wave  that  can  be  seen  just  leaving  the  region  at  the  bottom  of  the 
mesh.  This  situation  is  very  nearly  the  same  as  that  shown  in  Figure  53 
for  the  net  bypass  flow  case.  The  recirculation  eddy  is  stronger  here 
because  the  gap  flow  is  downward. 
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Figure  B-5.  See  caption  for  Figure  B-l.  Final  flow  pattern  is  a  simple  flow 
away  from  the  shell  and  downward  from  the  gap.  There  is  no  flow  from  the 
end  of  the  gap  in  this  frame  of  reference  but  the  scaling  of  vector 
lengths  accentuates  the  small  gap  flow.  The  flow  pattern  in  the  chamber  is 
very  similar  to  that  in  the  net  flow  case,  but  the  gap  direction  is 
reversed. 
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Figure  B-6t  Time  history  of  pressure  at  the  shell  corner  for  20  ram  diameter 
barrel  with  a  1%  gap  and  no  net  bypass  flow.  The  results  are  very  similar 
to  the  net  bypass  flow  case,  see  Figure  33,  because  the  pressure  level  is 
largely  set  by  the  prescribed  values  on  the  open  boundary. 
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Figure  B-7.  Time  history  of  temperature  at  the  shell  corner  for  a  20  mm  diam¬ 
eter  barrel  with  a  1%  gap  and  no  net  bypass  flow.  The  early  time  behavior 
is  similar  to  the  net  bypass  flow  case,  see  Figure  34,  but  after  about  0.7 
ms  this  case  produces  much  lower  temperatures.  This  is  caused  by  the 
cool,  expanded  gas  leaving  the  gap  and  mixing  with  the  wall  boundary 
layer . 


TEMPERATURE  CONTOURS  AT  TIME  =  .  I3496E-02 

Figure  B-8.  Temperature  contours  in  the  whole  field  visualization  for  a  20  mm 
diameter  barrel  with  1%  gap  and  no  net  bypass  flow.  The  temperatures  are 
considerably  less  than  in  the  bypass  case,  see  Figure  35,  because  there  is 
much  less  viscous  heating  in  the  gap. 
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TIME  =  . 22893E-C2  V  =  548.0 


Figure  B-9.  Whole  field  visualization  of  flow  pattern  in  40  mm  barrel  with  no 
gap.  See  Figure  B-l  for  description  of  distortion  introduced.  At  late 
times  flow  is  similar  to  20  mm  case;  see  Figure  4-9.  Early  times  showed 
the  previously  observed  rarefaction  motion  and  flow  reversal. 
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PRESSURE  CONTOURS  AT  TIME  ,22893E*02 

Figure  B-10.  See  caption  for  Figure  B-9.  Pressure  contours  at  late  time  are 
horizontal  and  nearly  constant. 
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Figure  B-ll.  Time  history  of  pressure  at  shell  corner  for  40  mm  diameter  bar¬ 
rel  with  no  gap.  Nearly  identical  to  20  mm  case,  see  Figure  10,  with 
characteristic  early  oscillations  due  to  rarefaction  propagation. 
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Figure  B-12.  Time  history  of  temperature  at  the  shell  corner  for  40  mm  diame¬ 
ter  barrel  with  no  gap.  Results  are  nearly  identical  with  20  mm  case. 

See  Figure  14.  Shell  mass  was  scaled  to  produce  the  same  accelerations. 
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Figure  B-13.  Whole  field  visualization  of  temperature  contours  for  40  ram,  no 
gap  case.  Strong  wall  boundary  layer  heating  is  similar  to  20  mm  case 
shown  for  corner  region  only  in  Figure  13. 
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Figure  B-14.  Whole  field  visualization  for  40  mm  diamete^  barrel  with  no  gap 
when  the  gas  kinematic  viscosity  is  assumed  to  be  10  cmz/sec.  See  Fig¬ 
ure  B-l  for  distortions  caused  by  this  plot.  Viscosity  is  about  an  order 
of  magnitude  greater  than  standard  case  but  flow  pattern  is  unchanged. 
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Figure  B-15.  See  caption  for  Figure  B-14.  Pressure  contours  at  late  time 

showing  the  horizontal  and  nearly  constant  values  characteristic  of  all  of 
the  no  gap  runs* 
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Figure  B-16.  Time  history  of  the  pressure  at  the  corner_of  t^e  shell  for  the 
40  mm  barrel,  no  gap  case  with  viscosity  equal  to  10  cmz/sec.  Behavior 
is  similar  to  other  no  gap  cases  except  that  rarefaction  oscillations  seem 
to  be  dampened . 
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Figure  B-17.  Time  history  of  temperature  at  the  shell  corner  for  the  40  mm, 
no  gap  case  with  viscosity  raised  to  10  cm  /sec.  Behavior  is  similar  to 
other  no  gap  cases. 
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Figure  B-18.  Whole  field  visualization  of  temperature  contours  at  late  time 
for  40  mm,  no  gap  case  with  viscosity  raised  to  10  cm  /sec.  Variation 
is  similar  to  that  seen  in  other  no  gap  cases. 
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Figure  B-19.  Whole  field  visualization  of  flow  pattern  for  40  mm  diameter 
barrel  with  a  1%  gap  when  no  net  outflow  is  permitted.  At  early  times  a 
rarefaction  region  develops  behind  accelerating  shell  and  gas  is  pulled  in 
from  gap.  Similar  to  20  ram  case;  Figure  B-l. 


TIME  =  . ! I705E-03  V  *  86.7 


Figure  B-20.  See  caption  for  Figure  B-19.  Recirculation  region  has  developed 
in  shell’s  wake.  Eddy  occupies  a  smaller  fraction  of  the  shell’s  base  than 
in  the  20  mm  case;  see  Figure  B-2.  Absolute  eddy  size  may  be  nearly  con¬ 
stant. 
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TIME  =  ♦ I 9383E-03  V  =  142.! 


Figure  B-21.  See  caption  for  Figure  B-19.  Recirculation  associated  with 
rarefaction  wave  is  very  similar  to  the  20  mm  case  shown  in  Figure  B-4. 
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Figure  B-22.  See  caption  for  Figure  B-19.  Final  flow  pattern  is  similar  to 
20  mm  case  in  the  chamber,  see  Figure  B-5,  but  differs  in  the  gap.  The 
wall  boundary  layer  is  unable  to  span  the  gap  so  countercurrent  flow  has 
developed  reminiscent  of  the  net  outflow  cases.  There  is  flow  splitting 
at  the  corner  in  this  case. 
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Figure  B-23.  Time  history  of  pressure  at  the  corner  of  the  shell  for  a  40  mm 
diameter  barrel  with  a  1%  gap.  Behavior  similar  to  20  mm  case  with  and 
without  net  bypass  flow. 
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Figure  B-24.  See  caption  for  Figure  B-19.  Pressure  contours  at  late  time. 
Pressure  in  chamber  exceeds  that  in  the  gap  in  order  to  accelerate  flow 
with  the  shell  and  against  the  wall  boundary  layer. 
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Figure  B-25.  Time  history  of  the  gas  temperature  at  the  corner  of  the  shell 
for  a  40  mm  diameter  barrel  with  a  1%  gap.  Large  oscillations  are  due  to 
vortex  shedding  at  the  corner  and  the  countercurrent  flow  into  the  gap 
which  did  not  existd  in  the  20  mm  case. 


TEMPERATURE  CONTOURS  AT  TIME  =  .229I7E-02 

Figure  B-26.  See  caption  for  Figure  B-19.  Temperature  contours  at  late  time 
show  values  comparable  to  the  20  mm  case  with  net  flow  bypass,  Figure  35, 
and  much  higher  than  the  20  mm  no  bypass  case,  Figure  B-8.  This  is  due  to 
the  countercurrent  gap  flow  which  brings  gas  into  the  gap  region. 
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Figure  B-27.  Whole  field  visualization  of  the  flow  pattern  for  a  40  mm  diame¬ 
ter  barrel  with  a  1%  gap  using  the  turbulence  model  of  Pope  and  Whitelaw. 
To  obtain  calculations  that  were  stable  the  turbulent  variables  had  to  be 
overwritten  whenever  they  exceeded  boundary  layer  values.  At  early  times 
the  flow  is  similar  to  the  standard  case;  Figure  B-19. 

TIME  =  .  I  I677E-03  V  =  86.5 


Figure  B-28.  See  caption  for  Figure  B-27.  Recirculation  eddy  in  shell's  wake 
is  similar  to  the  laminar  result  in  Figure  B-20. 
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Figure  B  29,  See  caption  for  Figure  B-27.  Reestablishment  of  recirculation 
caused  by  flow  reversal  of  a  rarefaction  wave  just  passing  out  of  the  sys¬ 
tem.  Compare  to  Figure  B-21. 
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Figure  B-30.  See  caption  for  Figure  B-27.  Boundary  layer  along  wall  much 
stronger  at  the  high  effective  viscosity  so  that  the  lower  gap  region  is 
pulled  downward  and  the  cross  flow  along  the  shell  bottom  is  deflected 
more  strongly. 
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Figure  B-31.  Time  history  of  pressure  at  the  corner  of  the  shell  for  a  40  mm 
diameter  barrel  with  a  1%  gap  using  the  Pope  and  Whitelaw  turbulence 
model.  Results  essentially  the  same  as  the  laminar  case;  see  Figure  B-23. 
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Figure  B-32.  See  caption  for  Figure  B-27.  Pressure  contours  for  whole  field 
showing  that  gap  pressure  is  less  than  laminar  case,  Figure  B-24,  in  order 
to  drive  gas  into  the  gap  against  the  -greater  viscous  forces. 
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Figure  B-33.  Time  history  of  gas  temperature  at  the  corner  of  the  shell  for  a 
40  mm  diameter  barrel  with  a  1%  gap  using  the  Pope  and  Whitelaw  turbulence 
model.  The  results  are  similar  to  the  laminar  case,  Figure  B-25,  although 
the  final  temperatures  are  higher  due  to  enhanced  viscous  heating. 


TEMPERATURE  CONTOURS  AT  TIME  -  .2284 !E *02 

Figure  B-34.  See  caption  for  Figure  B-27.  Whole  field  visualization  of  tem¬ 
perature  contours.  Values  are  about  150°C  higher  than  in  laminar  case, 
Figure  B-26,  because  of  enhanced  viscous  heating. 
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Figure  B-35.  Whole  field  visualization  of  flow  pattern  in  40  mm  barrel  with  a 
1%  gap  when  the  viscosity  has  been  raised  to  a  uniform  value  of  10 
cm2/sec.  See  caption  on  Figure  B-l  for  distortion  caused  by  this  kind  of 
plot*  Early  time  expansion  at  rear  of  shell  pulls  in  gasses  from  the  gap* 
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Figure  B-36*  See  caption  for  Figure  B-35.  Recirculation  eddy  at  rear  of 
shell  is  qualitatively  the  same  as  those  in  other  gap  cases. 
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Figure  B-37.  See  caption  for  Figure  B-35.  Reestablished  circulation  pattern 
following  rarefaction  wave.  Wall  boundary  layer  is  considerably  wider 
than  in  less  viscous  cases. 
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Figure  B-38.  See  caption  for  Figure  B-35.  Final  flow  pattern  is  similar  to 
other  no  gap  cases  in  the  chamber.  In  the  gap,  the  larger  viscosity  has 
strengthened  the  wall  boundary  layer  so  much  that  upflow  again  appears 
along  the  shell. 
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Figure  B-39.  Time  history  of  pressure  at  corner  of  shell  for  40  mm  diameter 
barrel  with  a  1%  gap  when  the  viscosity  has  been  raised  to  10  cm  /sec. 
The  behavior  is  similar  to  other  gap  cases. 


Figure  B-40.  See  caption  for  Figure  B-35. 
pressure  contours  showing  the  lowest  of 
low  pressures  develop  in  order  to  drive 


Whole  field  visualization  of 
all  pressures  in  the  gap.  These 
the  gas  into  the  gap. 
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Figure  B-41.  Time  history  of  the  gas  temperature  at  the  corner  of  the  shell 
f°£3a  ^  mm  diameter  barrel  with  a  1%  gap  when  the  viscosity  is  raised  to 
10  cm  /sec.  The  results  are  similar  to  other  gap  cases  but  temperatures 
are  generally  higher  due  to  enhanced  viscous  heating.  Irregularities  due 
to  flow  unsteadiness  are  now  more  evident. 


TEMPEPATURE  CONTOURS  AT  TIME  -  .229I0E-02 

Figure  B— 42.  See  caption  for  Figure  B- 35.  Whole  field  temperature  contours 
reveal  that  viscous  heating  has  raised  temperatures  400 °C  above  the  lami¬ 
nar  case. 
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Figure  B-43.  Expanded  computational  grid  for  calculation  at  standard  resolu¬ 
tion  (20  x  20  in  chamber).  Cells  have  been  shown  as  square;  correct 
aspect  ratio  is  shown  in  Figure  15. 
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Figure  B-44.  Whole  field  visualization  of  flow  pattern  for  a  40  mm  diameter 
barrel  with  a  1%  gap  and  a  viscosity  of  10  crn  /sec.  These  parameters 
are  the  standard  case  for  the  mesh  refinement  calculations.  Initial  flow 
separation  at  t  =  0.074  ms. 
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Figure  B-45.  See  caption  for  Figure  B-44.  Development  of  recirculation  zone, 
t  =  0.10  ms. 
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Figure  B-46.  See  caption  for  Figure  B-44.  Migration  of  eddy;  t  =  0.13  ms. 
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Figure  B-47.  See  caption  for  Figure  B-44.  Diffusion  of  eddy;  t  -  0,16  ms. 


TIME  =  . I 9380E -03  V  *  142.1 


Figure  B-48.  See  caption  for  Figure  B-44.  Flow  reversal  and  reestablishment 
of  recirculation  by  rarefaction  wave;  t  =0.19  ms. 
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See  caption  for  Figure  B-44.  Final  flow  pattern;  t  =  2.28  ms. 
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Time  history  of  pressure  at  cell  corner  for  calculation  at 
resolution. 
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PRESSURE  CONTOURS  AT  TIME  =  .22797E-02 

Figure  B-51.  See  caption  for  Figure  B-44.  Pressure  field  contours  for  stan¬ 
dard  resolution  run. 
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Figure  B-52.  Time  history  of  temperature  at  corner  of  shell  for  calculation 
at  standard  resolution. 
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TEMPERATURE  CONTOURS  AT  TIME  --  .22797E-02 

Figure  B-53.  See  caption  for  Figure  B-44.  Temperature  field  contours  for 
calculation  at  standard  resolution. 


Figure  B-54.  Expanded  computational  grid  for  making  the  twice  normal  resolu¬ 
tion  runs  (40  x  40).  Correct  aspect  ratio  of  cells  are  shown  in  Figure 
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Figure  B-55.  Whole  field  visualization  of  flow  pattern  for  (40  x  40)  resolu 
tion  calculation  of  40  mm  diameter  barrel  with  1%  gap.  Early  time  forma 
tion  of  rarefaction  region;  t  =  0.058  ms. 
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Figure  B-56.  See  caption  for  Figure  B-55.  Growth  of  recirculation  eddy; 
t  =  0. 12  ms. 
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Figure  B-57.  See  caption  for  Figure  B-55.  Growth  of  separation  eddy; 
t  =  0. 13  ms. 
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Figure  B-58.  See  caption  for  Figure  B-55.  Diffusion  of  recirculation  eddy 
t  =  0.16  ms. 
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Figure  B-59.  See  caption  for  Figure  B-55.  Reestablishment  of  recirculation 
eddy  due  to  rarefaction  wave;  t  =  0.19  ms. 


TIME  =  . 22863E -02  V  =  547.9 


Figure  B-60.  See  caption  for  Figure  B-55.  Final  flow  field;  t  =  2.29  ms. 
Note  that  growing  boundary  layers  on  the  shell  and  the  wall  lead  to 
countercurrent  flow  and  flow  splitting  at  the  corner. 
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Figure  B-61.  Time  history  of  pressure  at  the  corner  of  the  shell  showing  the 
same  behavior  as  in  the  lower  resolution  runs. 
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Figure  B-62.  Time  history  of  temperature  at  the  corner  of  the  shell  showing 
the  same  behavior  as  in  the  low  resolution  calculations. 
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Figure  B-63.  Computation  grid  for  highest  radial  resolution  calculation  (20  x 
80),  Grid  has  been  distorted  into  rectangular  configuration  for  the 
purposes  of  whole  field  visualization. 
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Figure  B-64.  Highest  resolution  calculation  of  40  mm  diameter  barrel  with  1% 
gap.  Early  time  expansion  region  behind  shell;  t  =  0.061  ms. 
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Figure  B-65.  See  caption  for  Figure  B-64.  Dissipation  of  recirculation  eddy 
t  =  0,12  ms • 
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Figure  B-66.  See  caption  for  Figure  B-64.  Dissipation  and  migration  of 
recirculation  eddy;  t  =0.15  ms. 
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Figure  B-67.  See  caption  for  Figure  B-64.  Reestablishment  of  the  eddy  by  the 
rarefaction  wave;  t  =  0.198  ms. 


Figure  B-68.  See  caption  for  Figure  B-64.  Late  time  flow  field  pattern; 

t  =  2.29  ms.  Strong  countercurrent  flow  in  gap  similar  to  40  x  40  resolu¬ 
tion  calculation. 
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Figure  B-69.  Time  history  of  pressure  at  the  corner  of  the  shell  for  40  mm 
diameter  barrel,  1%  gap  with  high  (20  x  80)  resolution.  Results  identical 
to  lower  resolution  runs. 
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Figure  B-70.  Time  history  of  gas  temperature  at  the  corner  of  the  shell  for 
case  described  in  Figure  B-69.  Results  substantially  the  same  as  in 
standard  resolution  calculation. 
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TEMPERATURE  CONTOURS  AT  TIME  -  .22S07E-O2 


Figure  B-71.  See  caption  for  Figure  B-64.  Temperature  contours  at  late  time. 
Values  agree  with  standard  resolution  case  to  within  about  20°C  except 
right  at  the  gap  exit.  See  Figure  R-53. 
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APPENDIX  C 


CONSIDERATIONS  ON  THE  OUTFLOW  BOUNDARY  CONDITIONS 
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Excerpt  from  letter  dated  27  February  1981 
to  Dr.  Csaba  Zoltani  from  Dr.  Stuhmiller 


It  was  a  pleasure  to  meet  with  you  and  the  BRL  staff  earlier  this  month  to 
discuss  both  the  flow  calculations  done  under  the  present  contract  and  the 
future  directions  of  BRL  in  numerical  simulations.  My  only  regret  is  that  the 
time  was  so  short  that  we  were  just  starting  to  get  to  the  nuances  of  the 
physics  and  the  numerics.  I  hope  that  over  the  next  several  months  we  can  rec¬ 
tify  that  situation. 

As  you  recall,  we  had  a  lengthy  discussion  on  the  flow  in  the  annular  gap.  In 
our  calculations,  we  assumed  that  the  average  gas  velocity  out  of  the  gap  was 
equal  to  the  shell  velocity.  You  expressed  your  concern  that  there  should  be  a 
net  positive  flow  rate  relative  to  the  shell.  When  we  considered  the  zero 
velocity  wall  boundary  condition,  then  we  find  a  region  where  the  velocity  is 
considerably  less  than  the  shell  so  that  the  net  flow  is  not  as  definite. 


The  curvature  of  the  profile  in  laminar  flow  is  governed  by  the  pressure 
gradient  and  can  take  a  maximum  value  up  to  the  sonic  speed.  It  was  decided  to 
arbitrarily  assume  a  flow  rate  that  was  10%  higher  than  the  shell  speed. 

Later  that  week,  Henry  Evans  was  contacted  about  using  an  integral  formula  for 
the  flow  rate  derived  from  the  steady  state,  one-dimensional  form  of  the  axial 
momentum  equation.  Taken  literally  this  formula  produces  supersonic  velocities 
(because  the  overall  pressure  drop  is  so  large  that  the  differential  form  of 
the  equations  is  incorrect).  More  importantly,  it  is  not  proper  to  apply  con¬ 
ditions  which  are  approximate  integrals  of  the  solution  -  the  partial  differ¬ 
ential  equatiions  must  have  local  boundary  conditions. 

Reopening  this  issue  led  me  to  look  into  this  compressible  channel  flow  more 
closely.  The  problem  can  only  be  resolved  with  a  two-dimensional  calculation 
with  shock  capture  if  it  exists,  of  the  channel  flow  and  the  widening  ahead  of 
the  shell,  but  this  goes  beyond  our  present  scope.  Instead,  I  considered  the 
allowed  flow  rates  for  laminar  and  turbulent  channel  flow  if  the  velocity  pro¬ 
file  maximum  was  sonic.  r 
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The  first  question  to  be  answered  is  whether  the  flow  is  laminar  or  turbulent. 
The  nominal  conditions  considered  are 

gap  half-width,  h  =  10  ^  m 

_o 

gap  length,  L  =  4  x  10  m 

o 

density,  p  =  200  kgm/m 

shell  velocity,  V  <  400  m/sec 

gas  viscosity,  y  «  8  x  10”"*  kgm/sec-m 


The  Reynolds  number  based  on  the  half  width  is 


Re 


»  io5 


which  suggests  that  the  flow  is  turbulent,  if  turbulence  has  time  to  develop. 
From  stability  theory,  disturbances  grow  as  e^ct  with  time  or  e^cx^  with  dis- 
tance.  For  plane  Poiseuille  flow  (according  to  C.  C.  Lin,  Theory  of  Hydrody¬ 
namic  Stability,  Cambridge,  p,  30)  the  maximum  growth  rate  occurs  at  Re  *  6  x 
10  for  a  wavelength  k  =  0.75/h  and  has  a  magnitude  c/U  »  7.65  x  10”^.  There¬ 
fore  disturbances  in  the  inlet  of  the  gap  will  grow  by  an  amount 


ekcL/V  _  e2.3 


by  the  end  of  the  channel.  This  growth  is  far  less  than  that  required  by  em- 

Q 

pirical  transition  theory,  the  so-called  e  rule,  so  despite  the  high  Reynolds 
number  the  flow  is  probably  unsteady  and  laminar,  not  turbulent.  These  con¬ 
siderations  are  the  same  ones  that  cause  me  to  doubt  if  true  turbulence  exists 
in  the  internal  ballistic  flows,  hence  if  conventional  turbulence  models  are 
appropriate. 

The  next  question  is  whether  laminar  flow  will  have  a  chance  to  establish 
itself,  that  is,  if  inlet  conditions  will  be  smoothed  out  at  the  parabolic 
profile  formed.  The  characteristic  time  behavior  for  the  slowest  viscous  mode 
is  e~Vt/ph2>  or  as  a  Spatial  variation  e"^x/pUh^#  Again,  over  the  length  of 
the  gap  the  inlet  profile  will  have  relaxed  by  an  amount  equal  to 


e-pL/pUh2  =  e-0.004 

Therefore,  by  viscous  effects  alone  the  inlet  condition  will  be  hardly  changed 
along  the  gap. 
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The  final,  simple  possibility  is  that  the  wall  boundary  layer  that  developed 
in  the  chamber  before  the  gap  triggers  plane  Couette  flow,  that  is,  a  linear 
variation  from  the  shell  velocity  to  zero  at  the  wall.  This  solution  is  dis¬ 
cussed  in  Liepmann  and  Roshko,  Elements  of  Gas  Dynamics,  Wiley,  pp.  311-328. 
Flow  separation  at  the  corner  is  probably  too  great  a  phenomena  to  ignore.  A 
schematic  diagram  of  the  corner  and  gap  flow  in  the  accelerating  coordinate 
system  is  as  follows 


77777777777^7777  ////// 


t~7~7~7  rrwv  -}  7-7-7  tt— 7-7— 7-7 


What  about  the  question  of  flow  rate?  We've  determined  that  the  flow  is  prob¬ 
ably  neither  fully  developed  laminar,  turbulent,  or  Blasius,  but  what  would  we 
get  if  it  were?  From  Schlichting,  Boundary  Layer  Theory,  McGraw-Hill,  p.  505, 
we  find  the  average  flow  velocity  (in  a  pipe)  is  related  to  the  maximum  veloc¬ 
ity  (which  we  assume  is  the  sonic  speed,  a)  by  the  formula 


u  /  u 
av  max 


(n  +  l)(2n  +  1) 


so  that  the  following  table  can  be  generated 


Re 


n 


max 


6.0 

6.6 

7.0 


0.791 

0.807 

0.817 


For  fully  developed  laminar  flow  the  profile  is 


u(y)  =  A  +  By  +  Cy^ 
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subject  to  the  boundary  conditions  that 


u(-h)  =  0 
u(+h)  =  V 


umax 


a 


which  after  some  simple,  but  lengthy  algebra  leads  to 


u  /u 
av  max 


where  M  is  the  shell  Mach  number,  V/a.  An  abbreviated  table  of  values  gener¬ 
ated  by  this  formula  is 


M 

uav^umax 

0.0 

0.667 

0.2 

0.698 

0.4 

0.725 

0.6 

0.744 

0.8 

0.749 

1.0 

0.667 

Finally,  plane  Couette  flow  would  have  a  profile 


u(y)  =  V  (l  +  y/h) 


with  the  corresponding  flow  rate  relation 


u  /u 
av'  max 


In  our  calculations  already  presented  we  used  the  condition 

uav^umax  ~  M 

The  various  flow  rates  are  displayed  on  the  enclosed  plot.  The  correct  flow 
rate  lies  somewhere  between  the  Couette  flow  that  would  exist  in  the  absence 
of  a  pressure  gradient  and  the  fully  developed  flow  which  we  know  do  not  have 
time  to  form.  The  assumption  used  in  the  present  calculation  is  in  that  broad 
band. 
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